Enhanced tunneling magnetoresistance in Fe|ZnSe double junctions 
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We calculate the tunneling magnetoresistance (TMR) of Fe|ZnSe|Fe|ZnSe|Fe (001) double mag- 
netic tunnel junctions as a function of the in-between Fe layer's thickness, and compare these 
results with those of Fe|ZnSe|Fe simple junctions. The electronic band structures are modeled by 
a parametrized tight-binding Hamiltonian fitted to ab initio calculations, and the conductance is 
calculated within the Landauer formalism expressed in terms of Green's functions. We find that 
the conductances for each spin channel and the TMR strongly depend on the in-between Fe layer's 
thickness, and that in some cases they are enhanced with respect to simple junctions, in qualitative 
agreement with recent experimental studies performed on similar systems. By using a 2D dou- 
ble junction as a simplified system, we show that the conductance enhancement can be explained 
in terms of the junctions energy spectrum. These results are relevant for spintronics because they 
demonstrate that the TMR in double junctions can be tuned and enhanced by varying the in-between 
metallic layer's thickness. 

PACS numbers: 85.75.-d, 72.25.Mk, 73.40.Rw, 73.23.Ad 
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A magnetic tunnel junction (MTJ) consists of two 
ferromagnetic electrodes separated by a thin non- 
conducting barrier. It is experimentally observed that 
the conductance of a MTJ depends on the relative ori- 
entation of the electrodes' magnetization, and because of 
this, during the last years a lot of attention has been paid 
to the investigation of MTJs as promising candidates for 
application in spintronic devices, such as read heads and 
magnetic random access memories (for reviews, see [1] 
and the references therein). 

One of the challenges, that has to be overcome for 
practical applications, is to reach higher values of the 
tunneling magnetoresistance ratio (TMR), defined as 
TMR= [(Tp-T A p)/Tp] x 100%, where T P and Tap are 
the conductances measured for the parallel (P) and an- 
tiparallel (AP) magnetization of the electrodes. Several 
possibilities are now being considered: to use highly po- 
larized materials (half-metals) or diluted magnetic semi- 
conductors as parts of MTJs, to produce junctions with 
almost perfect interfaces, and to use double magnetic 
tunnel junctions (DMTJs), in which metallic layers are 
inserted inside the semiconductor barrier of a MTJ. In 
this work we explore the latter alternative, and focus 
our attention on the dependence of the TMR on the in- 
between metallic layer's thickness. 

Since X. Zhang et al [2] suggested to use DMTJs, sev- 
eral groups [3-6] have theoretically shown that DMTJs 
exhibit richer spin-dependent transport properties than 
MTJs and that the TMR can be higher than that 
of MTJs, but only very recently could these DMTJs 
be fabricated [7,8]. T. Nozaki et al [7] have re- 
cently measured the tunnel magnetoresistance of epitax- 
ial Fe|MgO|Fe|MgO|Fe (001) DMTJs at room tempera- 
ture, and found an enhancement of the TMR with respect 
to MTJs (53 % for DMTJs versus 44 % for MTJs at low 



bias), indicating that DMTJs may present an advantage 
over simple junctions for their use in spintronics. 

As far as we know, up to now the only theoret- 
ical studies of DMTJs with magnetic layers in be- 
tween the semiconductor were made within the free elec- 
tron model (that cannot reproduce the decay rates in- 
side the semiconductor of evanescent states with dif- 
ferent symmetry), and using rectangular potential pro- 
files [2,3,5,14]. Moreover, these studies analyzed the de- 
pendence of TMR on the applied bias voltage and not 
on the in-between metallic layer's thickness, as we do 
in this work. For this reason, in this paper transport 
through Fe(oo)|ZnSe(a)|Fe(oo) (001) MTJs and through 
FeM|ZnSe(&)|Fe(c)|ZnSe(&)|Fe(oo) (001) DMTJs is the- 
oretically investigated using a realistic tight-binding 
(TB) Hamiltonian to obtain the electronic structure of 
the junctions. Fe(oo) are semi-infinite electrodes, and a, 
b and c denote thicknesses. The systems studied are epi- 
taxial, and we restrict to zero temperature, infinitesimal 
bias voltage and elastic transport. We choose Fe|ZnSe 
because it can be grown epitaxially and there is very lit- 
tle interdiffusion at the interfaces, thus producing crys- 
talline junctions in which there are no magnetically dead 
Fe layers [9,10]. Moreover, in contrast to what happens 
in FejMgO based junctions, there is no oxidation of the 
interfacial Fe layers, which is known to be detrimental to 
TMR [11]. To obtain a clearer insight into the physics 
involved in transport through double tunnel junctions, 
we also calculate the conductance through a simplified 
two-dimensional tunnel junction (2DDJ). 

The conductances are calculated from the active re- 
gion's Green's function G a s = [IE 

where 1 stands for the unit matrix, Hg is the Hamiltonian 
corresponding to the active region, are the self- 

energies describing the interaction of the active region 
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with the left (L) or right (R) electrodes (a corresponds 
to the majority or minority spin channels), and 'active re- 
gion' stands for whatever is sandwiched by the electrodes. 
For DMTJs, the active region consists of an 'in-between 
metal region' (IBMR) sandwiched by two identical 'semi- 
conductor regions' (SCR), while for MTJs the active re- 
gion is simply the SCR. The energy E is actually Ep+ir], 
Ep being the Fermi level of the system, and we take 
T] — > + . The self-energies are given by Yf L = H^ LS g1HLs 
and SJj, = HftggftHfts, where Hls and Hrs describe 
the coupling of the active region with the electrodes, 
and g a L j R are the surface Green's functions for each elec- 
trode. These surface Green's functions are calculated 
using a semi-analytical method [12] and are exact within 
our TB approximation. The transmission probability T a 
is given by [13] T a (k /h E F ) = Tr [AJGgA^G^] where 
^l/r = *(^L/fl — ^L/fl)' wnuc the conductance is given 
by 

" k n 

where Nk is the total number of wave vectors parallel to 
the interface that we consider (in our case 5000 is enough 
to achieve convergence in T). 

We start our discussion with the 2DDJs case, which are 
of the type M(oo)|S|M|S|M(oo), where M(oo) are semi- 
infinite paramagnetic metallic electrodes, S is a semicon- 
ductor and M is a metal (the same as the electrodes). 
The metal and semiconductor have the same structure, a 
square Bravais lattice with two atoms per unit cell, and 
are periodic in the direction perpendicular to the trans- 
port direction. The 2DDJs electronic structure is mod- 
eled by a 2nd nearest neighbors TB Hamiltonian with 
one s orbital per atom. The TB parameters are cho- 
sen to make Ep fall in the middle of the semiconductor's 
band gap (of 0.5 eV). The SC and IBM regions are varied 
between 3.2 A and 32 A . 

It is found that for certain thicknesses of the IBMR the 
conductance presents peaks in which it is enhanced by 1 
to 4 orders of magnitude, as can be seen in Fig. 1 for a 
2DDJ with a SCR of 12.8 A . This effect can be explained 
in terms of the active region's density of states (DOS), 
obtained from its Green's function Gs- When the con- 
ductance is enhanced, partial density of states (PDOS) 
calculations indicate that there exist states at Ep ex- 
tended throughout the whole junction, so in that situ- 
ations transport occurs through resonant states. When 
this happens, the 2DDJs conductances are higher than 
the corresponding ones of simple 2D junctions. These re- 
sults are consistent with those of Z. Zheng and coworkers 
for a DMT J with a non-magnetic in-between metal [14]. 

Fig. 2 shows the maximum ratio between the conduc- 
tance of 2DDJs and 2D simple junctions, as a function 
of the SCR thickness. The maximum attainable ratio is 
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FIG. 1: Conductance as a function of in-between metal thick- 
ness for 2DDJs with a semiconductor region of 12.8 A . 

of 146 % and occurs for a SCR thickness of 9.6 A and 
an IBMR thickness of 19.2 A . For thinner SCRs the 
ratio is nearly constant and roughly 140 %, but beyond 
12.8 A the enhancement effect is lost. Having mentioned 
the main results for the 2DDJs, we go on to discuss the 
details for the three-dimensional Fe|ZnSe DMTJs. 
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FIG. 2: Maximum ratio between the conductances of 2DDJs 
and simple 2D junctions as a function of the semiconductor 
region thickness. 

Fig. 3 shows schematically the structure of simple and 
double junctions, which are periodic in the x-y plane, and 
the different magnetic configurations considered, paral- 
lel (P) and antiparallel (AP). Since the coercive field of 
the electrode and the in-between Fe layers is different, 
the magnetic configurations shown are experimentally at- 
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tainable [7]. 
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FIG. 3: Schematic structure of simple and double junctions, 
and magnetic configurations considered: parallel (P) and an- 
tiparallel (AP). The junctions are periodic in the x-y plane, 
and the electrodes are semi-infinite. The arrows indicate the 
magnetization direction of the metallic regions. 

Fig. 4 shows the structure of a simple Fe|ZnSe junction 
with a SCR of 5.67 A , along the z direction (which is the 
direction of transport). The BCC Fe lattice parameter is 
2.87 A , and that of zincblende ZnSe is 5.67 A . 
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FIG. 4: Interfacial structure of a simple junction of BCC Fe 
electrodes and zincblende ZnSe semiconducting spacer along 
the (001) direction, which is the transport direction. In the 
example shown, the semiconductor's thickness is of 5.67 A . 

The electronic structure of the junctions is modeled by 
a parametrized 2nd nearest neighbors spd TB Hamilto- 
nian fitted to ab initio calculations [15,16], in which the 
hoppings between the Fe atoms and the (Zn,Se) atoms 
are calculated using Shiba's rules and Andersen's scaling 
law [17]. The Fe d bands are spin split by fiJdd, where 
ji = 2.2 fiB is the experimental magnetic moment of Fe 
and Jdd = 1-16 eV is the exchange integral between d 
orbitals (/xg is Bohr's magneton). With these values for 
fi and Jdd, the Fe d bands spin spitting is very well re- 
produced [15]. The ZnSe band structure is rigidly shifted 
to make the iron Fermi energy fall 1 eV above the ZnSe 
valence band and 1.1 eV below the conduction band, as 
indicated by photocmission experiments [9]. 

For simple junctions, we find that the conductances de- 
cay almost exponentially with semiconductor thickness, 
and that the TMR increases and is always positive (or di- 
rect), reaching a value of 90 % for a semiconductor thick- 
ness of 34 A . Our results arc in very good agreement with 
the ab initio results of MacLaren and coworkers [18]. 



For double junctions, we vary the SCR thickness be- 
tween 5.67 A and 28.35 A , and the IBMR thickness be- 
tween 2.87 A and 22.96 A . We find that the TMR and 
conductances strongly depend on the in-between metal 
thickness, and that for certain thickness combinations of 
the SC and IBM regions they can be higher than those 
corresponding to a MTJ, in agreement with the results 
of L. Shcng and coworkers [3]. The maximum ratio of 
DMTJs to simple MTJs conductances obtained is of 322 
%, and occurs for the P majority channel corresponding 
to SC regions of 22.7 A and an IBM region of 21 A . 
This large conductance ratio, which is pointing toward 
the existence of resonant states (confirmed by our DOS 
calculations), does not mean that the DMTJs TMR is 
going to be much larger than the MTJs one, although 
in general it is. In this particular case, the TMR value 
is of 97.9 % (the corresponding MTJ's value is 63.8 %), 
but in other cases the TMR values are greatly enhanced 
even in the absence of resonances. We find that the TMR 
enhancement can be a result of: (i) a drop in the con- 
ductance of some spin channels, while the conductances 
of other channels remain of the same order of magnitude 
as those in MTJs, or (ii) an increase in the conductance 
of one particular spin channel due to resonant tunnel- 
ing. Both effects are produced by a change in the active 
region's DOS near Ep, induced by the presence of the 
in-between Fe layers. 

As an example of resonance conductance enhancement, 
we show in Fig. 5 the active region's total DOS at Ep, as 
a function of the IBMR thickness and for a DMTJ with 
a SCR of 22.7 A . An increase in one order of magni- 
tude appears at an IBMR of 2.87 A for the P majority 
channel and for the AP minority channel, and a smaller 
increase appears at an IBMR of 8.6 A for the AP major- 
ity channel, while for the other cases the DOS is almost 
constant. These peaks coincide with a conductance en- 
hancement in these three channels, as it can be seen in 
the lower panel of Fig. 5, indicating that the origin of 
the conductance enhancement is the same as in 2DDJs, 
namely resonant tunneling. 

To visualize the interplay among the conductance val- 
ues of the different channels and configurations, Fig. 6 
shows the conductances and TMR values for a given 
DMTJ with a SCR of 17 A and those of the correspond- 
ing simple MTJ, as a function of the IBMR thickness. It 
is seen that, already for 6 A of Fe, the TMR is 1.5 times 
higher than that of a simple MTJ, although the conduc- 
tances are, in general, a little bit smaller. This is similar 
to what happens in all the cases studied. It is noticeable 
that for very thin Fe layers the TMR obtained for this 
DMTJ is negative, and that for IBMR thicknesses in the 
range 12-23 A the TMR is almost constant. This also 
happens for the other SCR thicknesses studied, and it 
is different to the damped oscillatory behavior that it is 
obtained using rectangular potential profiles and a non- 
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FIG. 5: Active region's density of states and conductances at 
Ef as a function of in-between metal thickness, for a semi- 
conductor region of 22.7 A . 
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FIG. 6: Conductance and TMR values of a Fe|ZnSe double 
magnetic junction with two semiconductor regions of 17 A , 
and those of the corresponding simple junction, as a function 
of the in-between metal thickness. 

For a given SC region thickness, we look for the max- 
imum attainable TMR value by sweeping over IBMR's 
thicknesses. We find that for ZnSe regions with thick- 
nesses below 20 A , the DMTJs' TMR values are 3 times 
higher than those of a simple junction, while the conduc- 
tances of some spin channels remain of the same order 
of magnitude. Beyond this thickness, the DMTJs TMR 
can be 50 times higher but negative (inverse TMR), al- 
though in this case the conductances of all spin channels 



are 4 to 6 orders of magnitude smaller, and thus very 
hard to measure. There is one particular case in which 
this docs not happen. For a DMTJ with a SCR of 22.7 
A and an IB MR of 8.6 A , we obtain a drop in the con- 
ductances of the P and AP minority channels and the P 
majority channel, and an enhancement of 175 % in the 
AP majority channel with respect to the corresponding 
MTJ, which results in a negative TMR enhancement by 
a factor of -40. 

In summary, we have investigated Fe|ZnSe double mag- 
netic tunnel junctions within a realistic Hamiltonian 
model and found that the TMR values can be much 
higher than those of simple junctions. We should mention 
that temperature effects, interfacial roughness, and the 
presence of defects in the DMTJs active region may de- 
crease the TMR values obtained in our calculations, but 
we believe that our results remain qualitatively valid. We 
conclude that the thickness of the in-between Fe layers 
in FejZnSe DMTJs is an interesting degree of freedom, 
which may make it possible to tunc and enhance the 
TMR of these systems, making them suitable for building 
future spintronic devices. To improve our understanding 
of these scarcely studied double junctions, it is highly de- 
sirable the experimental measurement of the TMR as a 
function of the in-between metallic layers thickness. 
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